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Abstract 

This paper studies the problem of robust downlink beamforming design in a multiuser Multi- 
Input Single-Output (MISO) Cognitive Radio Network (CR-Net) in which multiple Primary Users 
(PUs) coexist with multiple Secondary Users (SUs). Unlike conventional designs in CR-Nets, in this 
paper it is assumed that the Channel State Information (CSI) for all relevant channels is imperfectly 
known, and the imperfectness of the CSI is modeled using an Euclidean ball-shaped uncertainty 
set. Our design objective is to minimize the transmit power of the SU-Transmitter (SU-Tx) while 
simultaneously targeting a lower bound on the received Signal-to-Interference-plus-Noise-Ratio 
(SINR) for the SU's, and imposing an upper limit on the Interference-Power (IP) at the PUs. 
The design parameters at the SU-Tx are the beamforming weights, i.e. the precoder matrix. The 
proposed methodology is based on a worst case design scenario through which the performance 
metrics of the design are immune to variations in the channels. We propose three approaches based 
on convex programming for which efficient numerical solutions exist. Finally, simulation results are 
provided to validate the robustness of the proposed methods. 

Index Terms 

Robust beamforming, cognitive radio network, multi-user MISO communication, worst case 
design, imperfect CSI 
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I. Introduction 

A Cognitive Radio Network (CR-Net) Q, is an intelligent solution to the spectrum 
scarcity problem. In a CR-Net, the Secondary Users (SUs) are allowed to operate within the 
service range of the Primary Users (PUs), though the PUs have higher priority in utilizing 
the spectrum. There are two types of CR-Nets: opportunistic CR-Nets for which the SUs 
sense the spectrum and try to utilize the unused channels when they are not occupied by 
PUs; and concurrent CR-Nets in which SUs are allowed to use the spectrum even when PUs 
are active, provided that the amount of interference power to each PU is kept below a certain 
threshold J2|. Hereafter, we will focus on concurrent CR-Nets. 

Fig. [T]-a illustrates the downlink scenario of multiuser multiple-input single-output (MISO) 
CR-Net with K SUs coexist with L PUs. There are two sets of relevant channel state 
information (CSI) which play important roles in the system design: one set describes the 
channels between SU-Transmitter (SU-Tx) and SU-Receivers (SU-Rx's) while the other set 
describes the channels between SU-Tx and PU-Receivers (PU-Rx's). For simplicity, we term 
the first set of CSI as SU-link CSI and the second set as PU-link CSI. When PUs are inactive, 
the system becomes conventional multiuser MISO system, and SU-link CSI is needed for 
transmission design. This knowledge is usually acquired through transmitting pilot symbols 
from SU-Tx to SU-RXs, and feeding back the estimated CSI from SU-Rxs to SU-Tx. In 
practice, however, because of the time variant nature of wireless channels, it is not possible 
to acquire the CSI perfectly, either due to channel estimation error and/or feedback error. 
On the other hand, when PUs are active, PU-link CSI is further needed at SU-Tx for the 
purpose of controlling interferences at the PU-Rx's. This CSI knowledge has to be acquired 
by SU-Tx through environmental learning 0, which again will introduce errors in CSI. In 
this paper we consider the transmit design for a multiuser MISO CR-Net with uncertain CSI 
in both SU-link and PU-link. 

Previously in conventional radio network design, ad-hoc methods, such as diagonal loading 
JH, were exploited in the design procedure of robust beamforming systems. Quite recently 
these designs are based on well-reputed mathematical methodologies, such as the systematical 
worst case designs fl5)-|IE]|. These methods deal with a Minimum Variance Distortionless 
Response (MVDR) problem in the signal processing domain and show that the problem 
may be recast as a Second-Order Cone Program (SOCP) [9]. Also, it was shown that this 
worst case design scenario is equivalent to an adaptive diagonal loading [Q. One of the first 
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worst case designs was published by Bengtsson and Otterstten lfT0l -|[TT T| . They showed that 
the robust maximization of SINR would lead to a Semi-Definite Program (SDP) [9], after 
a simple Semidefinite Relaxation (SDR). Sharma, et al., |fT2j developed a model to cover 
the Positive Semi-definiteness (PSD) of the channel covariance matrix. They proposed two 
SDPs, a conventional SDP and a SDP based on an iterative algorithm. Also Mutapcic et 
al, lfT3l proposed a new tractable method to solve the robust downlink beamforming. Their 
method is based on the cutting set algorithm which is also an iterative method. Also |fl4 l |- lfT6l 
are targeting the robust design of a beamforming system using the worst case scenarios for 
Quality-of-Service (QoS) constraints. 

Quite a few works are published on the robust design for CR-Nets HH, H9l, [|20]| and EQ. 
Zhang et al. lfT8l . [fT9l have studied such a CR-Net from an information theoretic perspective. 
The CR-Net considered in JUS, 030 consists of one PU-Rx and one SU-Rx, and the SU-link 
CSI is assumed to be perfectly known, but the PU-link CSI has uncertainty. A duality theory 
was developed to cope with the CSI imperfectness. Additionally, the authors proposed an 
analytic solution for this case. Also, Zhi et al. GUI designed a robust beamformer for a CR- 
Net, where the system setup is the same as in lfT9ll . however there may be some uncertainty 
in both the channel covariance matrix as well as the antenna manifold. Finally, Cumanan 
et al. I12T1 considered a CR-Net having multiple PUs and only one SU. In this work, both 
channels are assumed to be imperfect. They also used the worst case design method to come 
up with a convex problem that can be solved efficiently. 

In this paper, we consider a downlink system of a CR-Net with multiple PUs and multiple 
SUs whose relevant CSI is imperfectly known. The imperfectness of the CSI is modeled using 
an Euclidean ball. Our design objective is to minimize the transmit power of the SU-Tx while 
simultaneously targeting a lower bound on the received Signal-to-Interference-plus-Noise- 
Ratio (SINR) for the SUs, and imposing an upper limit on the Interference-Power (IP) at the 
PUs. The design parameters at the SU-Tx are the beamforming weights, i.e. the precoder 
matrix. The proposed methodology is based on a worst case design scenario through which 
the performance metrics of the design are immune to variations in the channels. We propose 
three approaches based on convex programming for which efficient numerical solutions exist. 
In the first approach, the worst case SINR is derived through using loose upper and lower 
bounds on the terms appearing in the numerator and denominator of the SINR. Working 
in this line, SDP is developed which provides us the robust beamforming coefficients. In 
the second approach, the minimum SINR is found through minimizing its numerator while 
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maximizing its denominator. Different from the first method, we chose exact upper and lower 
bounds on the previously mentioned terms. This approach does not lead to a SDR but the 
resulting problem is still convex and may be solved efficiently. Finally, in our third approach, 
we find the exact minimum of SINR directly, and this method is also a general convex 
optimization problem. 

The rest of the paper is organized as follows. In Section HI] the model of our studied system 
is described. In Section HH] the robust design of a multiuser MISO CR-Net with multiple 
SUs and multiple PUs is considered. In Section [IV] we show that the resulting optimization 
problem, using loose upper and lower bounds, is a SDR In Sections |V] and |VT] we propose 
two more general problem formulations based on stricter bound and exact bound on the 
minimum value of SINR, respectively. In Section IVIII the simulation results that demonstrate 
the robustness of the proposed schemes are presented. Finally, in Section IVIIII we conclude 
the paper. 

Notations: Matrices and vectors are typefaced using slanted bold uppercase and lowercase 
letters, respectively. Conjugate and conjugate transpose of the matrix A are denoted as the At 
and A*, respectively. The trace of a matrix is annotated using Tr [■]. Positive semi-definiteness 
of the matrix A is depicted using A y 0. The symbol "=" means "defined as". C mxn is 
used to describe the complex space of m x n matrices. A Zero-Mean Circularly Symmetric 
Complex Gaussian (ZMCSCG) random variable with the variance of a 2 is denoted using 
CjV(0, a 2 ). For a vector like x, ||x|| is the Euclidean norm while the norm of a matrix like 
|| A || is the Frobenius norm. To show the differentiation of a function, /, with respect to 
some of its parameters, a, V a /(-) is used. Finally, mathematical expectation is described as 
E{-}- 

II. System Model 

Fig. Q] shows the downlink scenario of a multiuser MISO CR-Net coexisting with a Primary- 
Radio Network (PR-Net) having L PUs each equipped with a single antenna. The SU-Tx 
equipped with N antennas transmits independent symbols, s k , to K different single antenna 
SUs, {sk £ C}k=i- It is assumed that the transmitted symbols are all Gaussian, zero-mean 
and independent, i.e., s k ~ C7V(0, 1). Each symbol is precoded by a weight vector, {w fc G 
C Nxl }% = i, resulting in a vector signal, {s k = w fc Sfc}^ 1? for each one. It is known that 
s fc ~ CAf(Q, Qfe), where Q fe is the covariance matrix of s fc ; Q k = E{s fc s|,} = w&wj. y 
and is the zero vector. The channel from SU-Tx to each SU-Rx is determined using a 
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complex- valued vector, {h^ £ C Nxl }f =1 which is not perfectly known and there is some 
kind of uncertainty in channel gains. This uncertainty is described using an uncertainty set, 
Tik which is defined as an Euclidean ball 

W fc = {h|||h-h fc || <5 k }. (1) 

In this definition, the ball is centered around the actual value of the channel vector, h^, and 
the radius of the ball is determined by 5k which is a positive constant. Using this notion, the 
channel is modeled as 

h fc = h fe + a fc) (2) 

where is a norm-bounded uncertainty vector, ||afc|| < 5k- 

The SU-Tx combines the signals and transmits the combination, x, 

K 

x=^s Ji = Ws 1 (3) 

fe=i 

where s = [si, • • • ,sk] t £ C Axl contains the transmitted symbols and as we know, s ~ 
CSS (0, Ik), also W = [wi, • • • , w#] £ C NxK , is called the precoding matrix. For x we 
know that x ~ CSS(0, Q), where Q = E{xx^} = WW^ y 0. The design objective is to 
determine this precoding matrix W based on certain criteria that we will discuss in the next 
few paragraphs. 

The channel from the SU-Tx to a PU-Rx is also defined using a complex valued vector, 
i.e., {g£ £ C Nxl }f =1 . Here it is assumed that the CSI for these users is also uncertain. We use 
the same notation to describe the uncertainty for these channels. The uncertainty is defined 
using a set, Qi, which is 

Si = {s\\\s - Sell < Ve}- ( 4 ) 

Equivalently, we may write 

g e = g e + b e , (5) 

where h e is a norm-bonded uncertain vector, ||b^|| < rjt and ge is the actual value of the 
channel. 

III. Problem Formulation 
For the system depicted in Fig{T] the total transmitted power, TxP, is given by 

TxP^E{||x|| 2 } 

= EH w *ll 2 - (6) 

fe=l 
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The received signal at the A;th SU is 

y k = h{.x + n k 



K 



(7) 



i=i 

i^k 



The right-hand side of CD) has three terms. The first term is the received signal of the intended 
message, while the second and the third terms show the interference from other messages and 
noise, which is white and Gaussian, i.e. n k ~ £/V(0, cr 2 ), respectively. The average received 
power for /cth SU, S k , from the intended message is 

S k ±E{\hlw k s k \ 2 } 

= iwih fe | 2 . 

Similarly, it is easy to show that the received interference power, IP^, is given by 



(8) 



Eh 



k W i S i 



K 



E Iwjh, 



(9) 



i=i 

i^k 



Using ® and ©, the SINR of fcth SU-Rx, SINR fc , is given by 



SINR fc 



Iwth 



k"*\ 



i^k 



(10) 



In robust design problems relating to SINR, expressions with the form of |w£hfc| 2 are 
frequently used. We can write 

|w£h fe | 2 = w£(h fc + a fe )(h fc + & k )^ k 



= wl(H k + A k )w k , 
where H fc = h k h\ is a constant matrix and A k is given by 

A fe = h fc a{. + & k h{ + a fe a{.. 



(11) 



(12) 



Note that A k is a norm bounded matrix, || Afc|| < e k . It is straightforward to find the following 
relation: 



e k > II A A 



|h fe aj. + & k h{ + a*aj| 



< l|hjfeal| 



< llh* 



|a fc hl 



l a fe a ll 



^ + 25 fc ||h fe | 



(13) 
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Using (fl"3l) it is possible to choose e k = 5l + 2d>fc||hfc||. We may use the identity of x^Ax 



Tr 



Axx f 



, to further simplify this expression, which gives 

|w£h fc | 2 = Tr[(H fe + A fe )w fe wi] . (14) 

We also adopt the notation of = w^w]! in our design formulation. With this, we find 
that 

|wl.h fc | 2 = Tr[(H fc + A fe )w fc ]. (15) 

It is noted that, from now on, similar expressions will appropriately be used for the other 
terms of SINR fe . 

Also, the received signal at the Ith PU is 

z t = gjx + Vi 

K 

= J2ge w k s k + V£, (16) 

k=l 

and the interference power, IPi, to this PU-Rx would be 

K 
k=l 

Again using similar formulation as |w|.h fc | 2 , we get 

l w Ig^| 2 = W I( G ^ + A £ )w fc 

= Tr [(G, + Ae) W k ] , (18) 

where is a constant matrix, G^ = g^gj and is the norm bounded uncertainty matrix, 
IIA^H < Similarly we know that ^ = rj% + 2^||g^||. 

Our design objective is to minimize the transmitted power, TxP, while guaranteeing that the 
SINR at SU-Rx for all the channel realizations is higher than the QoS -constrained threshold, 
{SINRjt > 7i};[!Li> an d simultaneously IP at PU-PvX is less than the PR-net-imposed threshold, 
{IP£ < Ki}f =1 , respectively. Mathematically, this problem can be described as 

Minimize TxP (19) 

{w fe }f =1 

Subject to SINR fc > k = 1, • • • , K 

Vh k eH k 

TP t <K e , t = I,-- - ,L. 

The above problem is a problem with an infinite number of constraints. To deal with such 
a problem one well-known method is to find the minimum and maximum values of SINR fc 



August 5, 2009 



DRAFT 



!£,£.£. transactions on Wireless Communications (SuBmittea) 



and IPf, respectively, related to those realizations of the channels which are claimed as the 
"worst ones". The worst channel realizations for SINR and IP would lead to the minimum 
and maximum value of SINR and IP, respectively. In that case, the problem will guarantee 
that the smallest possible SINR and largest possible IP also satisfy the constraints. Using this 
worst case design methodology, we could recast (IT9T ) to a simpler problem set as follows: 



Minimize TxP 

jfc=i 



{w fc }f 



Subject to min SINR fc > j k , k = 1, • • • , K 
max IP^ < Kt, i = 1, • • ■ , L. 
Adopting the previously mentioned notations, we rewrite it as 

K 



Minimize Tr [W* 



{w fc }*" 



fc=i k=l 



Tr 



Subject to min 

l|A fc ||<e fc a2+Ef=iTr 



(Hfc + W k 



(H k + A k ) W, 



> 7 fc , k = 1, 



K 

max Tr \(G e + At) W k ] < Kt, £ — 1, ■ ■ ■ ,L. 

A/l \<nt , , L ^ ' - 1 



K 



(20) 



(21a) 



(21b) 



\\<nt k=1 

In the next sections, we will solve the robust problem of (|2T|) and will show that this problem 
can be recast as a series of simple optimization problems. 



IV. Loosely Bounded Robust Solution 

In this section we will deal with the problem of (|21"1) . In ifTOl and QT|, it is suggested to 
minimize the SINR through minimizing the numerator while maximizing its denominator. 
So (I21al) is equivalent to 



K 



min Tr 

l|A fc ||<e* 



^H fc + A fe )Wfc - 7fe max Tr ( H * + A fe) w i > 7fc^n- 



(22) 



As it is known, this method is a conservative way to find the minimum of the SINR. 

A. Minimization of SINR 
To minimize the numerator, 



min Tr 

l|A fe ||<e fc 



(23) 
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we adopt a loose lower bound, proposed by ifTOl . |fTTfl. Using this lower bound, we have 



min Tr 

l|A fc [|<e fc 



H fc + A fc W 



Tr 



(H fc - e k l N ) W» 



and to maximize the denominator, the following term should be maximized 



max Tr 

i|A fc ||<e fc 



(H* + A k ) W, 



Using a similar approximation, we have 

max Tr \(U k + A k ) W { 



Tr 



H k + e k I N W 



Using these results, the problem of SINR minimization (12 1 at) is recast as 



Tr 



K 

(fi fc - ejjv) W fc ] - 7fc E Tr [(H* + e fc Ijv) Wi] > o*7*, k = 1, - ■ ■ ,K 



8=1 



(24) 



(25) 



(26) 



(27) 



and by regrouping the left hand side of this equation we find 



Tr 



Hi 



K 



V 



i=l , 

i+k J 



efcTr 



K 



w fc + 7fc ^w l 



i=l 



><7= 7fc , & = !,•••, if. (28) 



B. The Whole Conventional Program 

Using the same methodology as before, IP maximization (I21bl) leads to the following 
problem 



K 

X;Tr[(G* + &I*) W A 



fc=i 



(29) 



Then the whole conventional program targeting to solve the robust downlink optimization 
in MISO CR-Nets becomes 

K 

Minimize Tr [W fc l 



Subject to Tr 



V 



K \ 

i^k J 



efcTr 



K 



1=1 

i^k 



Tr 



K 

(G, + ^Iiv)E W A 



k=\ 



<K £ , £= 1, 



W k = Wl k = l,---,K, 



W fc 0, fc = 1, 



> ^llk, fc = 1, • " - , K 

(30a) 
(30b) 

(30c) 
(30d) 
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Please note the fact that the last two constraints are inherent in the structure of the problem 
formulation. Also note that to come up with a convex problem formulation, a non-convex 
constraint, rankjW/c} = 1, is eliminated 0, ifTOl . |[T2l|. This final form of the problem is an 
SDP and can be solved using efficient numerical methods |[22|. Finally it should be noted that 
unlike lfT2l the beamforming weights are not exactly the principal eigenvector of the matrix 
solution. To get the beamforming weights, the eigen decomposition of the W k is used. In 
this decomposition, W k may be decomposed to a series of rank one matrices, i.e., 

N 

w k = J2 X n,k e n , k e{ jk , (31) 

n=l 

where in this expansion, X n>k denotes the nth eigenvalue and e n k is its respective eigenvector. 
The solution matrix of itself is a rank one matrix, then all the eigenvalues are equal to 
zero except one, let's say Xpf,k- Therefore the above mentioned equation may be written as 

Wjfc = Ajv.fc &N,k e jv,fc 

'^N,k e N,k)(\f ^N,k e N,k)^ 

= w fe wj,, (32) 

where = J\~N~ k ejv,fc- The next two sections deal with the same problem but in different 
ways. It should be noted that the next two problems are not SDP and are generally convex 
problems, but, for these problems we use the same formula to acquire the beamforming 
weights from the solution matrix. 

V. Strictly Bounded Robust Solution 

In the previous section, the minimum of SINR was found using loose upper and lower 
bounds for its constituent terms. In this section, we try to minimize the SINR using the same 
method: we minimize the numerator and maximizing the denominator. But here, we try to 
find the exact maximum and the exact minimum for each term respectively: 

K 

min Tr \{H k + A k )W k ] - lk ]T max Tr [(H fc + A fc )W,j > lk a 2 n . (33) 

l|A fe ||<£ fc L J f— J ||A fc ]|<e fc L J 

i^k 

Our main tool, is the Lagrangian Multiplier method. 

'The principal eigenvector of a rank one matrix is the eigenvector corresponding to the only non-zero eigenvalue. 
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A. Minimization of SINR 

We start with the first minimization problem. 



Proposition 1: For the terms Tr (H. k + A^J W fe , using a norm-bounded variable A k , 



|AJ| < efc, the minimizer and maximizer would be 



A mm 



w 



and 



W fc | 
t 



w, 

a max , fe 



(34) 



(35) 



respectively. 

Proof: Please refer to Appendix |A] 
Using the above results, we have 



min Tr 

i|A fc ||<e fe 



[H k + A k )W k 



max Tr 

|A fc ||<e fc 



So we may rewrite (1331) as 



= Tr 
Tr 



IWJ 











) w fe 


= Tr 


H/VV;, 


-Cfc||W fc 



H ' + e 'M w ' 



Tr 



H fc W, 



Tr 



H fc W fc 



i=i 



Tr 



V 



e fc ||W fc ||- 7fc ^(Tr H fc W 4 j+e fc ||W. 

/ 

V 



i=l , 

5. 77ze Whole Program 

Similarly, the IP constraints may be written as: 

K K 

„ m F E Tr f(G* + At) W fc l = E ( Tr 

l|A f ||<C fc fc=1 J fc=l 



|w fe ||+ 7fc EH w *l 



eJlWill. (37) 



(38) 



> 7fc<- (39) 



i=l , 



G e W k 



+ 6l|W fc || <«< 



(40) 
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Finally, the whole program is 



K 



Minimize Y^ Tr [WJ 



k=l 







( 


K \ 








A" 






Subject to Tr 




w fc 








l|W fc | 














i=l 1 








i=l 







fc = 1, 



fe=i 



+ 6 w fc < 



W fc = Wl, k = l,---,K 



W k hO, fc = 1, 



K. 



(41a) 
(41b) 

(41c) 
(41d) 



Although this final problem is not an SDP, it is in fact convex, and this problem can be 
solved using standard numerical optimization packages, like CVX ll22l . 



VI. Exact Robust Solution 

In the last two sections we dealt with the problem of minimizing the SINR using a 
conservative method. In this section we find the exact solution instead. We start again with 
the problem of (|2T|) . but with a simple alteration. This problem is stated as: 

K 

Minimize Y^ Tr [WJ 



Subject to min 

l!A fe ||< £fc 



A' 



Tr 



H fc + A fc )W fc - 7 fcE Tr (H fc + A fe )W 



8=1 



) 



K 



\\^d<tt 



k=l 



max Tr + A ^) W fc < K t 



L. 



>° n lk, k = l,---,K; 

(42a) 
(42b) 



In the above problem, we try to minimize the SINR directly and without using conservative 
assumptions. First we have the following proposition: 
Proposition 2: The minimizer of (|42al) has the form of 



a mm 
^k 



-ejfc- 



||Wfc- 7 fcE^i W» 

i+k 



(43) 



Proof: Please refer to Appendix El 
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Using this proposition we find 



mm 

|A fc ||<e fe 



A' 



Tr 



H k + A k )W k 



V 



Tr 



K 



- 7fc E Tr 



1=1 

i^k 



'H fe + A t )W, 



H t (W t - 7 ^W, 



i=l 
i^k 



K 



- efc ||W fe - 7*X] W *il - a n7fc- 



8=1 

i^k 



Finally we come up with the final and general problem: 



K 



Minimize Tr [WJ 



Subject to Tr 



K 



H t (W t - 7 ^W. 



8=1 



K 



K 



fc=i 



G,W fc 

t 



(44) 



- e fc ||W fc - 7 fe ^ W^l > a^7 fc , k = l,---,K; 

8=1 

(45a) 

+ 6l|W fe ||) < Kt, 



W k = Wl k = l,---,K 



W k hO, k = l,---,K. 



(45b) 

(45c) 
(45d) 



The above problem is the most general form of the original problem. It should be noted that 
the beamforming weights are also the principal eigenvector of the solutions of this problem. 
Also it should be mentioned that the IP part of these two last problems, (l4Tj) and (145]) . are 
the same, i.e., these two problems have the same performance in IPs. 



VII. Simulation Results and Discussions 

To validate our developed methods, a set of simulations were conducted. It is assumed that 
the BS is equipped with a Uniform Linear Array (ULA) having 8 elements with a spacing 
of half wave length. A set of K = 3 SU-Rx's are served and the CR-Net should protect a 
set of L = 2 PU-Rx's. The SU-Rx's are located in the directions of Q x = 20°, 9 2 = 35° and 
# 3 = 50° relative to the antenna boreside, respectively. The PU-Rx's are also located at the 
directions of 0i = 80° and 2 = 85°, respectively. It is assumed that the change in Direction 
of Arrival (DoA) of input waves to the SU-Tx may be changed up to ±5° arbitrarily. The 
noise power is assumed to be a 2 n = 0.01, and constant for all of the users. Also, a constant 
SINR level of lOdB is targeted for all the SUs, while the constant interference threshold of 
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0.01 is used to protect the PUs. The channel model for PUs and SUs is assumed to be in 
line with a simple model of 

[h*(0 fc )]i = e^-V™^, i, k = 1, • • • , K, (46) 
[&(<fit)]i = e i *<> i - r )<*»te\ i^=i,...,L. (47) 

The uncertainty sets are characterized with = r\t = 0.05. We have used the CVX Software 
Package ||22]| to solve the proposed problems numerically. 

In Fig. 12, the array gains toward each user using different weight vectors are depicted. 
In this figure and the subsequent figures, LCBS, SCBS and ExCS denote Loosely Bounded 
Convex Solution, Strictly Bounded Convex Solution, and Exact Convex Solution, respectively. 
In this figure, the vertical solid lines show the DoA corresponding to different SU-Rx's, 
while the vertical dashed lines show the DoA of PU-Rx's. From this figure, it is clear that 
the proposed method can transmit the desired data to the SU-Rx's while it is protecting the 
PU-Rx's. It is also apparent that all these approaches produce similar results. 

In Fig. |3l we have plotted the histogram of normalized constraints for both SU-Rx's and 
PU-Rx's. The normalized constraints of SUs, , defined as SINRk/lk, is equivalent to 

1 1 K 

c£ = ~r~ w i H fc w fc — ^5]wjH fc Wi, 

where H fc = h fc hj., and the normalized PU constraint, C% , is defined as 

cf p) = - E gjH fc g,. 

K t k=l 

Unlike the normalized SINR constraints for a SU-Rx, when a normalized IP constraint is less 
than one, this constraint is considered to be satisfied. In Fig. [3]-a and Fig. [3]-c, the normalized 
SINR histograms for two different scenarios are depicted. In the first one, the uncertainty 
sets are chosen to be = & = 0.05 while for the second scenario, the uncertainty measures 
are four times more than the first one, i.e. e k = £ e = 0.20. It is apparent that in this scenario 
the gap between the ultimate value of normalized SINR, i.e. 1, and the actual values is wider 
than in the first scenario, having smaller uncertainty sets. In both cases, as expected, ExCS 
is outperforming the other two schemes due to its exact bounds on SINR. In Fig. [3]-b and 
Fig. [3]-d the normalized IP constraints for PU-Rx's are depicted. As can be seen, there is 
little difference between the proposed methods in terms of their IP constraints, because of 
the similar structure of these constraints. As it can be seen, the robust design is immune for 
the variation of channel, whereas the non-robust model fails in such situations. Also it is 
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apparent that for the robust case, the variation of normalized constraints is much less than 
the variation of normalized constraints for the non-robust case. It is also clear that SBCS 
and ExCS are more efficient in terms of handling the SINR. The ExCS model not only can 
satisfy all the constraints, but also is a parsimonious model in terms of SINR. This is because 
of the fact that this method uses exact minimum of SINR. Also it should be noted that the 
IP variations in SBCS and ExCS are the same. It is because of the identical form of the 
equation which describes the IP constraints in these two methods. Additionally, it should be 
mentioned that they are slightly better than the LCBS. 

Also in Fig. |4] we have plotted the normalized total transmit power versus the SINR 
thresholds for different amounts of allowed normalized IP. The normalized total transmit 
power is the ratio of total transmit power to the noise power and the normalized IP is defined 
using the same manner. Both quantities are dimension-less and for better clarity are displayed 
in dBs. As expected, ExCS is better than the other two methods. In Fig. Hl-a, it is clear that 
for a relative IP level of —4 dB, ExCS transmits the lowest amount of power while SBCS 
requires to transmit a modest amount of power relative to LBCS, and finally LBCS requires 
to transmit the largest amount of power. In this figure, it is also observed that for a relative IP 
level of dB, the proposed ExCS is the best scheme to use to transmit power. In Fig.|U-b, we 
have plotted the same graph but in the higher SINR values. In this range of SINR thresholds, 
all of the optimization problems with a relative IP level of —4 dB would be infeasible, so 
the graph is only provided for the relative IP level of dB. It is clear that in such scenarios, 
ExCS performs the best, although it should be noted that the performance of SBCS and 
ExCS are very close to each other. 

VIII. Conclusion 

The problem of robust downlink beamforming design in multiuser MISO cognitive radio 
networks is studied. Particularly, a set up of K SU-Rx's and L PU-Rx's, all equipped with 
a single antenna is considered, and the SU-Tx has N transmit antennas. It is assumed that 
the relevant CSI is not known perfectly for both sets of users. The uncertainty in the CSI is 
modeled using an Euclidean ball notation. Three different approaches, namely LBCS, SBCS 
and ExCS, are presented which can be implemented efficiently. The first solution is a SDP, 
while the later two solutions are the convex optimization problems. Various simulation results 
are presented to evaluate the robustness of proposed methods. 
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Appendix A 
Proof of Proposition 1 

The Lagrangian function, using an arbitrary positive multiplier, A > 0, is 

L(A k} A) = Tr [(H fe + A k )W k ] + A(|| A fc || 2 - e 2 ) 



Tr 



[fl k + A k )W k + A(Tr 



(48) 



By differentiating II23II of this Lagrangian function with respect to A£ and equating it to zero 
we will find the optimizer A k , 



V A .L(A Jk ,A) = Wi + AA fe = 0, 



which is annotated by A° k pt , 



A Opt 
^k 



A k 



(49) 



(50) 



To eliminate the role of arbitrary parameter of A, again, we differentiate the Lagrangian 
function with respect to this unknown parameter and then equate it to zero 



V A L(A fc ,A) = 0, 
to get the optimizer A, annotated as X opt , 

A ? 3 * = -|| Wi||. 

By combining these results, finally, we come up with 



* opt 

k 



w 



(51) 



(52) 



(53) 



To test if this solution belongs to a minimum, we should observe that the second derivative 
at the optimizer points should have a non-negative value: 



V 2 ^L(A° k pt ,X opt ) = X opt > 0. 



(54) 



To find the maximum of such a term, again, using a positive arbitrary Lagrangian multiplier, 
we build a Lagrangian function. 



L(A k , A) = Tr (H fc + A fc )W t -A(||A fc || 2 -e 



Tr 



fH fc + AOW,- 



A(Tr 



el)- 



By differentiating it with respect to and equating it to zero 

V A *£(A fc ,A) = WJ-AA fc = 0, 



(55) 



(56) 
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we will get 



(57) 



Again, by differentiating the Lagrangian function with respect to A and equating it to zero, 



we are able to get the optimizer. 



V A L(A fc ,A) = 0, 



A°f* = -II Wdl, 



A opt _ 
k ~ e k 



|WJ 



To prove if this solution belongs to a maximum, we should observe that: 

V 2 ^L{A° k pt , X opt ) = -\ opt < 0. 



(58) 

(59) 
(60) 

(61) 



Appendix B 
Proof of Proposition 2 



The Lagrangian multiplier is adopted again: 



K 



L(A fc ,A) = Tr (H fc + A fc )W fc - 7 fcE Tr (H fc + A fc )W t + A(Tr 



i=l 
i+k 



By differentiating this function and equating it with zero, 



K 

V A ;£( A fc , A) = Wj - lk E W} + AA fe = 0, 



e\). (62) 



(63) 



=i 

i^k 



we will come up with 



Af 



1 

A 



K V 



V 



i=i 
i^k j 



and to eliminate the A, 



we will get 



V A L(A fe ,A) = ||A fc ||-e fc = 0, 



K 



A^ = _||W fe - 7fe E W ill' 



(64) 



(65) 



(66) 



i=l 
i^k 



A opt _ 
k ~ ~ e k- 



W fe - 7fc Ef=lW i 

V i-+k ) 

l|w fc - 7fc El=iW l 

i+k 



(67) 



The second order differential test to prove that this solution belongs to a minimum, in this 
case, is also straight forward and is not included here. 
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Fig. 1. A Typical Multiuser MISO CR-Net System with Uncertain CSI 
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Fig. 2. Array Gain for Different Users 
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